Quantum transport through a Tonks- Girardeau gas 
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We investigate the propagation of spin impurity atoms through a strongly interacting one-dimensional Bose 
gas. The initially well localized impurities are accelerated by a constant force, very much analogous to electrons 
subject to a bias voltage, and propagate as a one-dimensional impurity spin wave packet. We follow the motion 
of the impurities in situ and characterize the interaction induced dynamics. We observe a very complex non- 
equilibrium dynamics, including the emergence of large density fluctuations in the remaining Bose gas, and 
multiple scattering events leading to dissipation of the impurity's motion. 
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The impetus for miniaturization has resulted in the cre- 
ation of nanostructures in which the motion of particles is 
purely one-dimensional. In these systems motional degrees of 
freedom can be excited only along one direction whereas in 
the two orthogonal directions the system occupies the quan- 
tum mechanical ground state. To reach the one-dimensional 
regime the chemical potential and the temperature need to be 
much smaller than the transverse level spacing. Interacting 
particles confined to a one-dimensional wave guide are funda- 
mentally governed by many-body quantum mechanics [ 1 ] . In 
non-equilibrium situations this gives rise to genuine quantum 
dynamics, examples of which have been seen in single mode 
nanowires j^l and in atom traps B H Si- In this paper, 
we study the non-equilibrium transport of single or few impu- 
rity particles through a one-dimensional, strongly interacting 
Bose gas. The impurities are accelerated by a constant force, 
very much analogous to electrons subject to a bias voltage, 
and undergo scattering with the atoms in the Tonks-Girardeau 
gas. 

An interacting one-dimensional Bose gas realizes a bosonic 
Luttinger liquid. Its many-body quantum state in the ho- 
mogeneous case is characterized by a single parameter 7 = 
■p^^tzt Ist]. Here m is the atomic mass, giu is the ID cou- 
pling constant, and nio is the ID density. For weak interac- 
tions (7 ^ 1) Bose-Einstein condensation and superfluidity 
are possible in harmonically confined ID systems. For strong 
interactions (7 » 1) the longitudinal motion of the particles 
is highly correlated. In this so-called Tonks-Girardeau regime 
the Bose gas "fermionizes", i.e. its iV-particle wave function 
can be related to that of a iV-particle spin-polarized Fermi gas 
The density of the Bose gas as well as the den- 
sity dependent correlation functions become Fermion-like and 
superfluidity vanishes Iil2l 131. 

Both the weakly jsl and the strongly interacting ifTil [Tsll 
regimes of one-dimensional Bose gases have been accessed 
a few years ago. This experimental realization of the one- 
dimensional Bose gas with (5-functional interactions has trig- 
gered significant research efforts both experimentally and the- 
oretically 1 16]. Of particular interest have been dynamical ex- 
periments 1 3, 4, 5, 6]. Elementary transport experiments have 
seen the suppression of dipole oscillations in a conjugated po- 
tential [5] and the absence of thermalization |6]. These exper- 



iments, however, have focused on global properties of the gas 
rather than using single impurity probes. Moreover, in con- 
trast to previous non-equilibrium experiments in one dimen- 
sion, we work with an open quantum system in which the im- 
purity atoms continuously gain kinetic energy and can transfer 
this energy into the trapped gas by collisions. On the theoret- 
ical side, single particle perturbations have been studied in a 
number of different regimes illiflllTlEolEllElEllil. 




FIG. 1: (Color online) a) An array of one-dimensional Bose gases in 
a two-dimensional blue-detuned optical lattice. Vertical confinement 
is provided by a harmonic magnetic potential B{z). b) Creating im- 
purities by a radio frequency pulse resonant at a specific magnetic 
field. The impurity atoms in the |_F = 1, = 0) state experience 
no vertical confining potential and are accelerated by gravity into the 
-z- direction. 



Our realization of strongly interacting one-dimensional 
gases is depicted in Figure[T^. We start from producing almost 
pure Bose-Einstein condensates of *^Rb of up to 1.5 x 10^ 
atoms in a magnetic trap in the hyperfine ground state \F = 
1, mi? = —1). The harmonic magnetic trap has the frequen- 
cies llIx.z ~ 2tt X 39 Hz and ujy — 2tt x 11 Hz. Strong con- 
finement into a one-dimensional geometry is achieved by adi- 
abatically loading the three-dimensional Bose-Einstein con- 
densate into an optical lattice. The optical lattice is formed 
by two retro-reflected laser beams of wavelength A = 764 nm 
arranged in the horizontal xy-plane. In this blue-detuned opti- 
cal lattice the atoms are trapped in the intensity minima of the 
interference pattern. Thus the vertical confinement is purely 
magnetic. At the position of the condensate, the standing 
wave laser fields overlap perpendicularly with orthogonal po- 
larizations and are focused to a circular waist (1/e^-radius) of 
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180 /im. The frequencies of the two beams are offset with re- 
spect to each other by 280 MHz. The optical potential depth 
U is proportional to the laser intensity and can be expressed 



in terms of the recoil energy Erec = 



2m>? 



Adiabatic load- 



ing into the ground state of the optical lattice was achieved by 
ramping up the laser intensity with a linear ramp of 150 ms du- 
ration. Typically, we confine approximately 50 atoms per one- 
dimensional tube in the center of the three-dimensional gas 
cloud. The interaction parameter 7 depends on the strength of 
optical lattice and ranges between 2 and 10. Here 7 represents 
the density-weighted average of across all tubes. The 

error in the determination of 7 is 15% dominated by uncer- 
tainties in atom number and calibration of the lattice depth. 

The hybrid magnetic/optical trapping potential provides us 
with the necessary means to create and detect the impurities 
with very good spatial resolution. We create spin impurities 
by using spatially resolved radio frequency manipulation. We 
apply a pulse (200 /is) of radio frequency resonant with the 
|F = l,mF = -1) -^\F = l,mF = 0) transition B . The 
spatial width of our addressing region is Fourier-limited by 
the duration of the pulse and is Az « 2.3/im (see FigurefTJ)). 
We drive approximately a 7r/2 pulse producing an impurity 
wave packet containing up to 3 atoms per one-dimensional 
tube. The atoms in the impurity state have zero magnetic mo- 
ment and are accelerated downwards by gravity. They start 
from zero center-of-mass velocity and reach several times the 
velocity of sound at the edge of the trapped cloud. During 
the propagation they are continuously confined to the one- 
dimensional waveguide in the radial direction because the ki- 
netic energy acquired is less than the radial level spacing. On 
their way downwards they strongly interact with the remain- 
ing atoms because the three-dimensional scattering lengths for 
collisions in the three possible combinations of states are ap- 
proximately equal to each other 

In order to accurately study the motion of the impurity 
atoms in our one-dimensional Bose gas we have performed 
a time resolved tomographic measurement of both the Bose 
gas and the impurity density distribution. At a variable time t 
after preparation of the impurity wave packet we measure the 
density distribution in each of the components separately in 
situ. We employ magnetic field (position) sensitive microwave 
transitions (pulse duration 200 /is) between the hyperfine 



ground states \F — l,Tnp = —1) \F — 2, m^? = 0) 
(trapped gas) and \F — l,mF = 0) ^ |F = 2, = 1} 
(impurity) and detect the atoms by absorption imaging on the 
\F = 2) ^ \F' = 3) transition of the D2 line. Figure|2]shows 
three snapshots of the time evolution of both the trapped com- 
ponent (blue) and the impurity (red). We observe that the im- 
purity atoms are initially well localized in a compact wave 
packet whose width is in agreement with the Fourier limit. 
As the impurity wave packet propagates it becomes wider and 
distorted indicating a strong dispersion and dissipation. For 
example, at (r =2000 /is, z=-10 /im) we observe a very dis- 
tinct asymmetric steepening of the propagating impurity wave 
packet which resembles the proposed shape of a supersonic 
shock wave in one dimension [26]. Moreover, the propaga- 
tion of the impurity leaves a strong imprint on the trapped 
component indicating a strong mutual interaction. We observe 
that the statistical noise of the data points in the sample after 
propagation of the impurity (r =1400 fis and t =2000 /is) is 
considerably larger than for the initial state. This could be at- 
tributed to high frequency oscillations of the atomic density at 
a length scale below our spatial resolution limit. 

We have analyzed the time evolution of the center-of-mass 
and the width of the impurity component (see Figure|3]l. Very 
clearly, the center-of-mass of the impurity does not follow 
a ballistic trajectory but its motion is hindered by the pres- 
ence of the interacting Bose gas. Theoretical investigations 
for impurities of equal mass exist only for the case of non- 
accelerated impurities JllllHi^ IH SlUlil and mostly 
for low momenta of the impurity. In our case, however, the 
impurity atoms are accelerated by gravity and reach the ve- 
locity of sound c already after traveling a distance shorter 
than the interparticle separation I/tiid- In the conceptually 
simplest case of binary collisions the probability for a mo- 
mentum change in a collision between two atoms in ID (a 
reflection of the impurity particle from the majority atom) is 
p = {1 + {aiomv/hfy^ |27]. Here am = 2?iV(™.9ic) 
is the one-dimensional scattering length. For the first colli- 
sion (i.e. after traveling a distance of order l/ni^i) we find 
p Ki 0.2 which becomes quadratically smaller for subsequent 
collisions as the velocity of the accelerated impurity particle 
increases. Consequently, in the majority of the binary colli- 
sions the impurity is transmitted. For binary collisions the en- 
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FIG. 3: (Color online) The circles show the measured center-of-mass 
position taken for 7 = 7. The error bars are the statistical error of the 
center of mass of the measured density distribution. The solid line is 
the prediction according to the model described in the text. The gray 
shaded area indicates the regime of uncertainty of 10% of uid given 
by our experimental parameters. The dashed curve indicates purely 
ballistic motion. The squares show the increase of the width of the 
impurity wave packet. The data point at 2 ms contains atoms which 
have already left the trapped gas which is not taken into account by 
the theory. 



ergy d issipation scales like Ebin{v) oc /{I + {mai^jv/h)'^) 
OalSTj]. Refining this result for a many -body system we adopt 
an approach based on Fermi's golden rule. The rate of scat- 
tering r{ki, kf) a single impurity particle from an initial mo- 
mentum state hi to a final momentum state kf can then be de- 
termined by the dynamic structure factor S{q, ui) of the gas: 
T{ki,kf) oc J dujS{q,uj) i28il . Here q — ki — kf is the mo- 
mentum change of the impurity, huj — e{ki) — e{kf) is its 
corresponding energy loss, and ujq is the excitation spectrum 
of the gas. We assume a free particle like dispersion relation 
for the impurity atoms (which is viable to first order in the 
impurity-gas interaction strength |22]). For a homogeneous 
weakly interacting system we have calculated the energy dis- 
sipation rate of the impurity using E{ki) oc J dqduoLoS{q,Lo) 

H to be Esf{v) = ^^St^v (1 - (f )') for v > c. In 
a homogeneous Tonks gas we find the energy dissipation 
Etg{v) = ^ma"^" ^ ~ —amv for u > vp ^ hnnio/m, 
i.e. a constant, density dependent force. In the limit of large 
velocities both results agree with the result for binary col- 
lisions. In figure [3] we plot the predicted trajectory from a 
numerical integration of the equation of motion taking into 
account a Thomas-Fermi profile of the density of the major- 
ity component as well as antibunching over a length scale 
of We use the predicted value of the force constant 

a{z) — with a{z = 0) = — 8m/s^ and find very 

good agreement with the experimental data without any ad- 
justable parameters. This theoretical model neglects the dis- 
tortion of the density of the trapped gas as well as coherent or 
multiple scattering processes and back action onto the impu- 
rity atoms. 

The impurity wave packet spreads considerably as it prop- 
agates. We have measured the root-mean-square (rms) radius 



w of the wave packet and plot its evolution in Figure [3] A 
function w{t) — ^/w^+vH^ is fitted to the data to extract 
the velocity v with which the wave packet spreads while it 
is fully inside the trapped gas. We find w = 1.5 ± 0.4mm/s 
which is slightly smaller than the uncertainty limited momen- 
tum distribution of the impurity wave packet of 2.0mm/s. 
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FIG. 4: (Color online) Time delay between the impurity preparation 
and the exit of the last impurity atoms from the trapped cloud for 
different lattice depths. The starting position is relative to the center 
of the trapped cloud. The solid lines show fits according to the model 
described in the text using a Thomas-Fermi profile for the density. 
The numerical prefactor agrees with the model to within a factor of 
3. 

We have investigated the number of collisions the impuri- 
ties are undergoing during the transport and the time lag of 
their motion. Using state-selective detection we determine 
how the number of impurity atoms inside the trapped com- 
ponent decays as a function of time. From these data we ex- 
trapolate the time delay after which all impurity atoms have 
left the trapped Bose gas for different initial positions and for 
different values of 7 (see Figure |4]l. In a simple model we 
assume that in every collision event the impurity motion is 
set back to zero velocity and afterwards the atoms are accel- 
erated again by gravity. On average the time between two 
collision events can be estimated tcou ~ \/2/ {niD{z)g) in 
which g is the gravitational acceleration. The delay time 
as compared to ballistic motion can be expressed as the to- 
tal number of inelastic collisions times the time delay accu- 
mulated per collision event (approximately tcoii) resulting in 
Jl°j,n,D{^)Tiz)t,ouiz? dz. Herer(z) = ^^^^^ is the 

collision rate for v « y^2g/niDiz) > c. This simple model 
explains the observed behavior well (see fits in Figure|4]i. The 
number of collisions is on the order of 5 < t delay /tcoii ^ 20, 
depending on the interaction strength and the starting position. 

Figure|5] shows absorption images of impurities propagated 
out of a Tonks gas for different values of 7 while the one- 
dimensional confinement was kept on. The image is recorded 
a fixed time after the preparation of the impurity wave packet. 
Impurity atoms leave the Tonks-Girardeau gas after a variable 
time depending on how many collisions they undergo. Af- 
ter the atoms have left the sample they simply undergo bal- 
listic motion. We observe a leading wave packet followed 
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FIG. 5: (Color online) Time of flight image of the impurity atoms 
and corresponding profiles of the optical column density. The im- 
ages are averaged over 5-10 repetitions of the experiment. Impuri- 
ties released from a lattice of potential depth 15 i5,.ec (7 = 3) (a, 
squares), 30£'rec(7 ~ 5) (b, open circles), and 45-Erec(7 ~ 7) (c, 
full circles) taken 18.7 ms after the radio frequency pulse. Both the 
number of atoms scattered out of the main peak and the length of the 
distribution of the scattered atoms increases with increasing 7. The 
curves in (d) have been offset vertically for clarity. The width of the 
atomic density distribution corresponds to the first Brillouin zone of 
the optical lattice. 



by a relatively long tail of atoms. From the position of the 
main peak we conclude that the atoms have not undergone 
momentum changing scattering while moving through the 
one-dimensional sample. The width of the main peak corre- 
sponds to an rms velocity spread of (1.9±0.1) mm/s, in agree- 
ment with the uncertainty limited width of the impurity wave 
packet. The atoms in the tail are scattered out of the main 
peak and, qualitatively, the tail represents the one-dimensional 
analog of the s-wave scattering spheres previously observed in 
three dimensions ll29ll but including multiple scattering events. 
The length of the tail is in good agreement with the results pre- 
sented in Figure |4] We observe a distinct minimum between 
the main peak and the tail (see Figure |5}l) for large values 
of 7. This minimum could reflect the increasing fermionic 
nature of the atoms in the gas. Since our impurity initially 
was a constituent of the Tonks-Girardeau gas, the two-body 
collision rate may still be related to the density-density corre- 
lation function g^'^\r) and becomes suppressed on the length 
scale of r K, l/niD as the gas fermionizes |fl2il . Then col- 
lisions between the impurities and the trapped atoms would 
be suppressed at the very beginning of the dynamics. Quan- 
titatively, the width of the minimum of scattered atoms in the 
time-of-flight image agrees with a suppression of collisions 
over the length scale « l/rtiu in situ. A comparison exper- 
iment performed with a weakly interacting three-dimensional 
Bose-Einstein condensate did not reveal a similar feature. 



In conclusion, we have studied quantum transport of 
spin impurity atoms through a strongly interacting one- 
dimensional Bose gas. For pseudospin-1/2 or spin-1 Bose 
gases as the majority component our tomographic detection 
technique could reveal fundamental properties of spin trans- 
port II19II and spin-charge separation [20.1 . 
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